Project1D340

library(igraph)

## 
## Attaching package: 'igraph'

## The following objects are masked from 'package:stats':
## 
##     decompose, spectrum

## The following object is masked from 'package:base':
## 
##     union

library(plotly)

## Loading required package: ggplot2

## 
## Attaching package: 'plotly'

## The following object is masked from 'package:ggplot2':
## 
##     last_plot

## The following object is masked from 'package:igraph':
## 
##     groups

## The following object is masked from 'package:stats':
## 
##     filter

## The following object is masked from 'package:graphics':
## 
##     layout

library(ggnetwork)

df <- read.table("~/Downloads/higgs-reply_network.edgelist")

ig <- graph_from_data_frame(df)

degs <- data.frame(ind=(ig %>% degree(mode="in") %>% log(base=10) ))
degs$outd <- ig %>% degree(mode="out")
in_hist <- ggplot(degs, aes(x=ind)) + geom_histogram()
ggplotly(in_hist, text=ind)

## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

## Warning: Removed 20332 rows containing non-finite values (`stat_bin()`).

ebs_df <- data.frame(abs=edge_betweenness(ig))

g <- ggplot(ebs_df) + geom_histogram(aes(x=log(abs,base=10)))

subg <- subgraph.edges(ig,1:500, delete.vertices = T)
tiny_g <- ggnetwork(subg, arrow.gap=0,layout=nicely())
tg <- ggplot(tiny_g, aes(x=x,y=y,xend=xend,yend=yend)) + geom_nodes()
tg <- tg + geom_edges(arrow=arrow(length=unit(4,"pt")))
tg